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ABSTRACT 

The effect of environment on galaxy formation poses one of the best constraints on the inter- 
play between mass assembly and star formation in galaxies. We present here a detailed study 
of the stellar populations of a volume-limited sample of early-type galaxies from the Sloan 
Digital Sky Survey, across a range of environments - defined as the mass of the host dark mat- 
ter halo, according to the groups catalogue of Yang et al. The stellar populations are explored 
through the SDSS spectra, via projection onto a set of two spectral vectors determined from 
Principal Component Analysis. This method has been found to highlight differences not seen 
when using standard, model-dependent comparisons of photo-spectroscopic data. We find the 
velocity dispersion of the galaxy to be the main driver behind the different star formation his- 
tories of early-type galaxies. However, environmental effects are seen to play a role (although 
minor). Our Principal Components allow us to distinguish between the effects of environment 
as a change in average age (mapping the time lapse of assembly) or the presence of recent star 
formation (reflecting environment-related interactions). Galaxies populating the lowest mass 
halos have stellar populations on average ~1 Gyr younger than the rest of the sample. The 
fraction of galaxies with small amounts of recent star formation is also seen to be truncated 
when occupying halos more massive than Mjj ^ 3 x 1O 13 M . The sample is split into satellite 
and central galaxies for a further analysis of environment. Small but measurable differences 
are found between these two subsamples. For an unbiased comparison, we have to restrict 
this analysis to a range of halo masses over which a significant number of central and satellite 
galaxies can be found. Over this mass range, satellites are younger than central galaxies of 
the same stellar mass. The younger satellite galaxies in Mjj ~6x 1O 12 M halos have stellar 
populations consistent with the central galaxies found in the lowest mass halos of our sample 
(i.e. Mh ~ 1O 12 M ). This result is indicative of galaxies in lower mass halos being accreted 
into larger halos. 

Key words: methods: statistical - galaxies: elliptical and lenticular, cD - galaxies: evolution 
- galaxies: formation - galaxies: halos. 



1 INTRODUCTION 

The origin and evolution of early-type galaxies is a long debated 
topic, its solution involving a large array of cosmological and as- 
trophysical processes. The current paradigm of galaxy formation is 
embedded in the ACDM cosmology, from which the structures in 
the Universe are built hierarchically. Under this framework, early- 
type galaxies are effectively a secondary stage in galaxy evolution. 

* E-mail: ferreras@star.ucl.ac.uk 



The first stage consists of the formation of rotationally-supported 
disk galaxies, built up through the accretion of gas and smaller sys- 
tems. Mergers subsequently operate, creating early-type galaxies. 
The masses and ages of the stellar populations in early-type galax- 
ies imply that these systems are the result of an intense starburst, 
followed by processes which quench star formation after which the 
galaxy evolves passively. 

One of the proposed mechanisms to stop star formation re- 
quires the removal of gas from the galaxy, involving either fast re- 
moval of cold gas (ram-pressure stripping), or a slower removal 
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of hot diffuse gas (strangulation). These mechanisms, however, do 
not result in a change in kinematic s and only minor changes in 
morphology (Weinmann et al. 2009) and so they can possibly ex - 
plain the increased fraction of SO galaxies (e.g. lDressleiJll98(ib . 
The major formation process of red s equence galaxies is t hought 
to operate through major mergers (e.g. lDe Lucia et al.l2 006) which 
result in the required structural and dynamical changes. Such 
processes have bee n known for a considerable time to produce 
spheroidal galaxies dToomre & Toomrdl 19721 : iBarnes & Hemquistl 
1996; iKhochfar & Burked 120031) , a s well as more detailed pho - 
tometric and dynamical properties dNaab & Burkertll200i 12006b . 
More recent results also suggest that minor mergers may play an 
increasingly important role in the build up and size evolution of 
massive ellipticals at relatively later times (e.g. Khochfar & Silk 
l2006l : lBezanson et al.ll2009l : IBernardi et alj [2009). The subsequent 
quenching of star formation and evolution of the galaxy onto the 
red sequence, requires the gas to be removed or heated to pre- 
vent the formation of new stars. Currently models invoke feed- 
back from active galactic nuclei (AGN), since this fits naturally into 
the merger scenario. Such interactions are expected to drive mate- 
rial to the centre of the galaxy through tid al torques, towards th e 
central supermassive black hole (SMBH) ( di Matteo e t al. 2007). 
The discovery of a correlation betwee n the mass of the SMBH 
and galaxy mass dGebhardt et alJ koOO). put significant weight be- 
hind the idea and provided scope for more comprehensive theo ries 
dHopkins et a"i]|2006l : iFaber et alj|2007l : ISomerville etai]|2008l) . 

This scenario naturally introduces an expectation of environ- 
mental dependence, since such formation process involves inter- 
actions with neighbouring galaxies and structures. A correlation 
with environment can arise in two forms. Firstly through the initial 
conditions, as these provide the impetus for the formation of the 
first galaxies so that objects in dense enviro nments will form ear- 
lier than in average or low density regions (Gottlobe r et al.ll200ll ; 
lBerlindetal]|2003b . Secondly, in higher density regions interac- 
tions, mergers, gas stripping, etc, are more likely to take place 
over the lifetime of a galaxy and so galaxies in these environ- 
ments will be pushed onto the red sequence at earlier times. Cer- 
tainly the fact that red early-type gala xies are preferentially found 
in higher density environmen ts (e.g. [Dressier 1980; Blant orTet al.l 
120051 ; Weinmann et al. 2006), suggests that environment plays an 
important role in their evolution. Therefore, looking at environ- 
mental differences in the stellar populations of early-type galaxies 
offers a method by which to constrain their formation. 

There has been much work on studies of differences in stel- 
lar populations through many different methods; the variations in 
the tight correlations followed by the early- type population such 
as the colour magnitude relation ( see e.g. Gallazzi et alj [2006b 
and the fundam ental plane (see e.g. iBernardi et alj 2003b , galaxy 
colour s (see e.g. lBlanton et alj|2005b , absorption line indices (see 
e.g. lNelanetalJl2005b and the parameters of population synthe - 
sis modelling (see e.g. IBernardi et ail 120061 ; iThomas et ai]|2005h . 
In all cases the effect of environment has been shown to be rela- 
tively weak, if observed at all. Thus we take a different approach, 
choosing a different methodology involving principal component 
analysis on spectral data to identify small d ifferences betwee n the 
stellar populat ions of early-type gal axies ^Rogers et alj|2007b . in a 
similar style to lFerreras et alj d2006h . over a range of environments. 

The environment is in most cases quantified through the pro- 
jected number density of galaxies, typically the distance to the n th 
nearest neighbour. However it has been argued that more physically 
motivated scales of environme nt are the mass and the virial radius 
of the host dark matter halo |Kauffmann et al. 2004; lYang et all 



120051 ; IWeinmann et al.ll2006l ; iBlanton & Berlin j|2007l) . Not only 
are environmental dependencies observed to act primari ly over dis- 
tance s comparable to the virial radius of such halos dGoto et all 
2003), but also the merger history o f the dark matter halo is de- 
termined mainly by its present mass ( Kauffmann et al. 200|). The 
mass of the host dark matter halo cannot be directly measured 
in most cas es but can be est imated through galaxy group cata- 
logues (e.g. lYangetaT]|2007b . Such catalogues also provide the 
halo-centric radius and can be used to easily separate the sam- 
ple into central and satellite populations. The other advantage in 
estimating the environment through halo mass is that it allows a 
direct comparison betwe e n observations and the o retical models. 
iDekel & Birnboiml j2006b.lSomerville et alj j2008b . lCattaneo et alj 
(2008) or lKhochfar & Ostriker! d2008b all make model predictions 



of the evolution and properties of gala xies in terms of the dark mat- 
ter halo mass. For ex ample Cattaneo et al. (2008), following th e 
work of iBirnboim & Dekell d2003h and lDeke l & Birnboiml d2006h . 
suggested that the downsizing observed in elliptical galaxies can 
be modelled by considering a critical mass halo above which gas 
cannot be accreted efficiently, being shock heated to the virial tem- 
perature, thus effectively shutting down star formation. 

This paper is structured as follows: we describe the sample of 
early-type galaxies used in this study as well as the details of the 
principal component analysis. We investigate the results of the PC A 
projections over a range of halo mass and velocity dispersion, we 
highlight the differences o bserved and investigate th em using the 
stellar population models of lBruzual & Charlotl|20 03). The sample 
is finally split into central and satellite galaxies, whose properties 
are compared. The satellite population is then used to determine a 
possible dependence on halo-centric radius. 



2 THE SAMPLE 

This work is based on th e large sample of early-type galax- 
ies of | Rogers et al] 12007). This sample is selected from the 
IBernardi et al. 1 2006b catalo gue, comp i led fr om the Sloan Dig- 
ital Sky Survey (SPSS lYorketall l2000b . Data Release 4 
dAdelman-McCarfhv et alj [2006). It is a volume-limited sample 
within z^0. 1 and M r ^ —21. A cut with respect to signal to noise 
ratio (S/N) was also imposed, rejecting those spectra with S/N^ 15. 
The final sample compri ses 7,134 early type galaxies. Here we ex- 
tend lRogers et alj d2007b to investigate the effect of environment in 
more detail, including the information of the host dark matter halo 

for each galaxy, a explained below. 

W e make use of the Galaxy Groups Catal ogue of Yang et all 
( 120071) , which is an improved application of the lVang et al.1 J2005) 
halo-based galaxy group finder to the New York University Value- 
Added Catalogue dBlanton et al.ll2005h, based on the Sloan Digital 
Sky Survey Data Release 4 dYork et a 1 2000). We cross correlate 
this catalogue with the original sample to find halo masses for all 
but 175 galaxies, leaving a total of 6,959 galaxies in the new sample 
used here. We have also removed a small (~ 150) set of galaxies 
with velocity dispersions below 125 km/s, as these galaxies only 
appear in the bin with the lowest halo masses, and have no counter- 
parts in more massive halos. Allowing them to be part of the lowest 
velocity dispersion bin would considerably reduce the average a in 
this bin, introducing a bias when compared to other hal o mass es. 



The galaxy group finder, described in I Yang et al.1 (2005) and 
lYang etalJd2007h . is an iterative process in which the membership 
of the groups and the relationships between the properties of the 
halo are refined at each step. Initially the group finder uses a friends 
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Figure 1. Our sample of SDSS early-type galaxies is shown with the two 
main parameters that are used to characterise the intrinsic (velocity disper- 
sion; a) and environmental (host halo mass; Mh) dependence of the under- 
lying stellar populations. In the crowded regions of the figure, contour lines 
track the density of galaxies in the plot. The grid corresponds to the binning 
applied throughout the paper. The large black dots and error bars give the 
average and RMS within each bin in halo mass. 

of friends algorithm with small linking lengths to identify the cen- 
tres of possible galaxy group candidates. In such groups the centre 
of the halo is given by the luminosity-weighted centre. The remain- 
ing isolated galaxies not associated with groups are also set as the 
potential centres of groups. The characteristic luminosity (Lig.5) 
of each candidate group is then estimated, where this is defined 
as the summation of the luminosity of all group members with 
01 M r — 5 log h ^ — 19.50- The L19.5 values are corrected for sur- 
vey completeness and the apparent magnitude limit of the survey at 
redshifts z^0.09. 

Using the characteristic luminosity, the mass of the group halo 
is estimated from a group mass-to-light ratio. The ratio used is set 
at a constant value across all groups for the first iteration but is 
subsequently refined to be group mass dependent. However, the re- 
sults are not particula rly sensitive to t he exact value of the ratio 
even if it is held fixed (Yang et al. 2005). The estimate of the group 
halo mass allows the derivation of other group halo properties such 
as the halo radius, within which the halo has an average density 
contrast of 180, and the virial radius, defined as the radius within 
which the average density is above a set value. Once the proper- 
ties of the h alo have been estimated, a NFW profile is used for the 
dark matter ( iNavarro. Frenk & Whitel ll997) to determine the three 
dimensional density contrast of the halo in redshift space. Further 
galaxies are subsequently assigned to the galaxy group candidates 
if they are within a certain distance of the centre. This process - 
which allows for the merging of two groups if all members sat- 
isfy the above criteria singularly - is repeated until the membership 

1 01 M r is the SDSS-r band magnitude for which K and E corrections are 
applied at redshift z=0. 1 



Figure 2. The projected radial distance from the luminosity-weighted centre 
of the group is shown with respect to velocity dispersion. The distance is 
scaled to the virial radius, as determined from the properties of the halo. 
The grid overlaying the sample shows the binning applied in the analysis. 
Black dots and error bars correspond to the average and RMS scatter of the 
sample within the bins in R pr oj ■ 

of each group remains constant. The final dark matter halo mass 
is then estimated from a linear relationshi p with respect to stellar 
mass, derived from semi-analytic models (Kang et al. 2005]). The 
galaxy group differentiates between central and satellite galaxies. 
Those galaxies which are the most massive galaxy of the group are 
defined as centrals, whereas the remaining galaxies are labelled as 
satellites. Low mass halos can consist of a single, central galaxy. 
Shown in figure [TJ is the entire sample as a function of host halo 
mass, Mh, and central velocity dispersion, a, of the individual 
galaxies. The black dots correspond to the average and RMS scatter 
of the sample, within the halo mass bins shown by the grid. 

We can see from the figure that within the lowest halo mass 
bin (Mh ~ 10 12 ft M©) there is a limitation imposed by the size 
of the halo on the highest sigma galaxy that it can contain. This is 
similar to that seen in the main catalogue in reference to the max- 
imum stellar mass (see e.g. lYang et al.ll2008T) . The difference here 
is the non-trival mapping of velocity dispersion to dynamical or 
stellar mass, which is dependent on other factors such as structure. 
However note that since the results presented in this paper are con- 
sistent across the range of velocity dispersions considered, the de- 
tails of such mapping would be unlikely to affect them. Also we are 
using velocity dispersion as a measure of the intrinsic properties of 
the galaxy, which is consistent with many stellar population s tudies 
in the literature (e.g. lBernardi et alj2003l : lGailazzi et al.ll2005h . 

The figure also reveals a shift towards decreasing velocity dis- 
persion at lower halo mass, whereby less massive halos on average 
contain smaller galaxies. This effect is illustrated by the solid cir- 
cles and error bars in figure [TJ which represent the average and root 
mean square deviation of the velocity dispersion within each bin of 
halo mass. 

This correlation b etween mass / luminosity an d environment 
in terms of halo mass jvan den Bosch et al]|2008j) but also seen 
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Figure 3. Left: Distribution in PCA components r) and £ of the upper and lower bins in velocity dispersion (bottom, labelled in km/s) and host halo mass 
(top, labelled in log(M/MQ)) of the entire sample in r\ and f. Right: distribution of the KS statistic (Dks) f° r me subsamples of the same bins in velocity 
dispersion (dashed) and halo mass (long-dashed) for r) (Bottom) and £ (Top). In order to quantify the statistical significance of the segregation, the result of 
drawing subsamples at random over all mass ranges and environments is shown (solid). 



with respect to galaxy density dBlanton et alJl2005l ; ICroton et al.l 
2005), means that it is important to take into consideration how 
such intrinsic properties change as a function of environment. It is 
well known that the s tellar populations of a galaxy are related to 
its velocity dispersio n ( Kauffman net alj2004| -| Thomas et al .120051 ; 
iGallazzi et a"fl2 005), therefore such a correlation can generate spu- 
rious environmental effects. Hence, in order to carry out a robust 
analysis, it is important to separate out the effects of environment, 
from those that are caused by this selection bias. 

To achieve this, we investigate the difference with respect to 
environment only at the same a range. Our sample is divided into 
sub-samples either with respect to halo mass, or velocity disper- 
sion. This classification is shown as a grid in figure Q] The choice 
of the grid size is motivated by the density of the underlying sam- 
ple, selecting larger bins at high mass. 

We also consider the variation in the properties of stellar pop- 
ulations with respect to the position of a (satellite) galaxy within 
its host halo. The dependence is measured as a function of the pro- 
jected distance from the luminosity weighted centre of the group, 
namely the projected halo-centric radius, R pro j- This will allow 
us to assess the effect of galaxy accretion in groups on the star 
formation history. We expect galaxies on the outskirts of a group 
to be newer members. Hence, these galaxies will be increasingly 
subject to the various effects of the group environment, such as 
'strangulation', 'ram-pressure stripping' a nd 'harassment' which 
act on the satellite galaxies of the groups iWeinman n et aljfeoOol ; 
Ivan den Bosch et al T l2008dlbT) . Shown in figure [2] is the satellite- 
only sample as a function of projected distance scaled by R v i r . 
Similar studie s of satellites have found a radial mass segregation 
within groups ( Ivan den Bosch et al . 2008a), with the least massive 
galaxies at the outskirts of the group. We find no such trend within 
our sample with respect to velocity dispersion. 



3 THE STELLAR POPULATIONS OF DIFFERENT 
HALO MASSES 

Regarding the effect of environment on the stellar populations of 
elliptical galaxies, it is important to notice that the reported dif- 
ferences have generally been small. For example, studies of the 
colour-magnitud e relation have only f o und limited and stati stically 
weak evidence ( iBernardi et"all 120031; IGallazzi et al.l [2 006). even 
more comprehensive analyses Jciemens eta l. 2006; B ernardi et al.l 
2006) find difference of the order ~ 1 Gyr. In this paper we op- 
timise the extraction of differences from spectroscopic data via 
principal component analysis, which has been shown to succeed 
in detecting differences within highly h omogeneous samples (e.g. 



in detecting ditterences withm highly h 
iFerreras et al.||2006l ; lRogers et al. 2007) 



3.1 PCA 

Principal Component Analysis (PCA) is a multivariate technique 
that reduces the dimensionality of a data set. In most cases this is 
just used as a data compression algorithm. However, previou s work 



impri 

dMadgwick et alj2003l ; lFerreras et alj2006l ; lRogers et alj2007l) has 
shown that vital information can be extracted from the projections 
on to the first few principal components. In this work the variables 
that describe the data set are the flux values at each wavelength 
i.e. the spectra. The task of PCA is to generate a set of basis vec- 
tors (the principal components) from the data set, such that one 
can rank these vectors with respect to the variance they capture. 
Hence, when obtaining the "coordinates" of each galaxy by pro- 
jecting their spectra on to the principal components, one can use 
just the very few coordinates that corr espond to the princi pal com- 
ponents with the highest variance. In Rog ers et alj J2007h we find 
that the first two components already hold valuable information re- 
garding the average age of the stellar populations and the presence 
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of recent star formation. We refer the reader to that paper for a de- 
tailed description of the methodology, although a brief summary is 
given below. 

The first principal component is found to be consistent with 
a typically old stellar population, showing a pronounced 4000A 
break, significant metal absorption lines and limited Balmer line 
strengths. The second component is much bluer, with absorption 
lines dominated by the Balmer series. A correlation between both 
principal components is forced by the orthogonality inherent to 
PCA. When projecting onto the observed spectra, the components 
give a relation that represents the mass-metallicity-age correla- 
tion o f early-type galaxies (e.g. lBernardi et all2003l ; lThomas et al.l 
2005). The extended scatter in the direction of the second compo- 
nent i.e. with respect to an excess of blue light, is created by recent 
star formation. This issue is confirmed in two ways: 1) through the 
application of a two component stellar population model, where 
we find galaxies with a higher projection onto the second prin- 
cipal component require a higher mass fraction in young stars 
jRogers et ai]|2007l . 120091) . 2) comparing the results of PCA with 
NUV photometry from GALEX: The NUV— r col our is highly sen- 
sitive to small amounts o f recent star formation ( Sch awinski et al.l 
l200dlKavirai et alj2007h . Within our sample, 'NUV bright' galax- 
ies (NUV— r$;4.9) present higher projections of the second princi- 
pal component, compared to the sample of quiescent, 'NUV faint' 
galaxies (NUV— r^5.9). This trend between PCA and the presence 
of young stars is optimised by rotating the projections on the two 
dimensional plane spanned by PC 1 and PC2, to give two new com- 
ponents: 77 and f . These two components can be defined as the dis- 
tance along the PC1-PC2 correlation (77) and the residual of the 
correlation (Q. Even though NUV is more sensitive to the presence 
of massive stars than optical spectra, this method is complementary 
to NUV studies. An advantage of using optical spectra over NUV 
photometry is that we can track the presence of recent star forma- 
tion for a greater length of time: NUV light decays very rapidly as 
the mo st massive sta r s die o ut. 

In lRogers et alj j2007) we compare the PCA projections with 
a number of star formation histories combined with population syn- 
thesis models and conclude that 77 is mainly related to the average 
age of the stellar populations, whereas C tra cks the presence of re- 
cent star formation (see also iRogersI 12009). This interpretation is 
carried over to this paper for the analysis of the effects of environ- 
ment. 



3.2 Stellar Mass vs Group Halo Mass 

In this section we look into the dependence of the stellar properties 
of elliptical galaxies on both the stellar mass and the mass of the 
halo which the galaxy occupies. We utilise the results from PCA 
described above, focusing on how the average properties of the 
galaxy (through rf) and the young stellar populations (through Q 
depend on galaxy mass or halo mass. While it is true that environ- 
ment plays a major role on the number de nsity of early- type galax- 
ies (i.e. the morphology-density relation, lDresslenll980l) . here we 
pose a different question: "at a fixed a, how different are the star 
formation histories of early-type galaxies with respect to environ- 
ment?". The left panel of figure [3] shows the distribution of the first 
and last bins in both velocity dispersion (bottom) and group halo 
mass (top). The comparison reveals that in terms of the mass of 
the galaxy (i.e. velocity dispersion), smaller galaxies have higher 
values of both 77 and ( indicating a younger age, most likely due 
to increased r ecent star formation , in agreement with the 'downsiz- 
ing' scenario dCowieetalJI 19961) . However the distributions of the 
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Figure 4. The average of the r\ and f components as a function of a within 
each group halo mass bin, as labelled (given in log(M/MQ )). The error bars 
correspond to the error on the mean. The small panels below the main plots 
show the standard deviation of r\ and £ in each bin. The thin black line is 
for reference only and is a least squares fit to the three most massive halos. 



two extreme bins regarding environment (top) show a considerably 
smaller divergence in terms of their 77 or C, distributions. Note that 
the preliminary results show that those galaxies in the lowest mass 
halos have slightly higher values of 77 and 

The use of a Kolmogorov-Smimov (KS) test confirms whether 
the galaxies in these subsamples originate from the same distribu- 
tion. This is quantified in terms of the Dks statistic (a measure of 
the maximum di fference of the cumulative distributions, see e.g. 
IPress et all 2009). The right panel of figure [3] shows the result of a 
Monte Carlo simulation, where we extract samples of 100 galax- 
ies from the upper and lower mass bins, and perform the KS test 
in each case. This exercise is repeated 10,000 times, and the his- 
togram of Dks is shown when extracting galaxies from the upper 
and lower bins in velocity dispersion (short dashed line), in host 
halo mass (long dashed line) of from random sampling of the com- 
plete sample (solid line). Velocity dispersion clearly plays the dom- 
inant role; its histograms in the 77 and £ components are far from the 
distribution of random sampling. However, environment - or halo 
mass - plays a more subtle role, although a non-negligible one, es- 
pecially with respect to the 77 component, which does show less of 
an overlap with the distribution for random sampling. Whether this 
effect is driven by the mass bias mentioned above or whether is it a 
genuine effect of environment is studied in this paper. 



3.3 A detailed look at the [77, Q distributions 

Many authors have studied the role of en vironment on the scal- 
ing relationships. T his includes colours (e.g. lWeinmann et al.l2006l: 
iKavirai et alj|2007 ), absorption line indices (e .g iKuntschner et all 
120011 : iNelan et alj|2005l : iBernardi et"ai1l2006l) as well as derived 
quantities such as luminosity- weighted age, metallicity and [a /Fe] 
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Figure 5. The fraction of galaxies with a £ value above 0.5, which is con- 
sistent with the presence of recent star formation. Each line corresponds to 
a range of host halo masses, as labelled (given in log(M/MQ)). The error 
bars are given by Poisson errors. 



(Thomas et al. 2005; Clemens et al . 2006). In a similar vein, we ex- 
tend this analysis to include the 77 and £ components. 

Shown in figure [4] are the average values of 77 and £ with re- 
spect to velocity dispersion, for a range of host halo masses. Points 
corresponding to the same host halo mass range are connected. The 
velocity dispersion value for each bin is given by the average. A ref- 
erence line is also plotted in the bottom panel. It is a least squares fit 
to the most massive halos and is shown in all subsequent figures of 
(77) for comparison. The error bars give the uncertainty on the mean 
values within each bin. The RMS deviation is shown separately at 
the base of each panel as cr(C) (top) and a{r\). 

A clear dependence is found with respect to velocity disper- 
sion, such that galaxies with lower values of a have higher compo- 
nents 77 and reveali ng the presence of yo unger stellar populations. 
Note that as shown in Roge rs et alj fc007t) . principal components 77 
and C are not dependent on the velocity dispersion, and so no cor- 
rection is needed (in contrast with analyses of equivalent widths). 

It is important to notice from figure [4] that the effect of envi- 
ronment is very subtle. In terms of (£) and in most cases (77), the 
trends of the different host halo mass ranges are indistinguishable 
from each other. However, galaxies in the least massive halos show 
a consistent and significantly higher value of (77), across all values 
of a (although the nature of the trend is the same for all halo mass 
bins). Given that the differences are small, the significance is deter- 
mined through a KS test for the 77 component, performed by com- 
paring galaxies in the lowest and the highest host halo bin, across 
all values of a. In the three central bins of velocity dispersion we 
find a probability higher than 99% that the samples are drawn from 
different distributions. At the highest and lowest bins of a the small 
number of galaxies prevent us from stating a similar result. 

Furthermore, figure [4] also indicates that (() does not differ- 
entiate with respect to halo mass. This comes as a surprise, as one 
might have expected that recent star formation is responsible for the 
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Figure 6. Change in the 77 component relative to the average relationship 
with a via SSP fitting (see text for details). A comparison of the lowest 
mass halo (dashed line) to the most massive halos (solid line) is shown 
relative to deviations from this relationship due to an increase/decrease in 
the average SSP age across cr (dotted lines). 



higher values of {77). To explore this issue in more detail, we use 
the conditional fraction of galaxies with a value of ( above a thresh- 
old at which one needs to invoke recent star formation. We choose 
£ ^0.5, since galaxies above this value both require significant 
fractions of young stars in a two-burst model an d also have NUV lu- 
mino sities consistent with recent star formation ( Roge rs et alj2007l 
120091) . 

In figure [5] we show this conditional fraction as a function of 
velocity dispersion for a range of host halo masses. Consistent with 
the simple analysis performed above, the main driving force behind 
the fraction of galaxies with recent star formation is galaxy mass 
(assuming a is representative of the mass of the galaxy). With re- 
spect to halo mass, the fraction of younger galaxies seem to split 
such that there is a consistent lower fraction for the most massive 
halo bins (black solid line) relative to the rest of the sample. The 
significance of this drop is not high in any one bin but the trend 
is notably consistent a cross all bins. This is a q ualitatively similar 
result to that found in ISchawinski et al I j2006f) using NUV data, 
although here we see that recent star formation persists in halos be- 
low Mh ^ 3 x 10 13 /i -1 Mq and is reduced in more massive halos. 
The possible mechanisms underlying this effect are discussed in the 
conclusions. 



3.3.1 Modelling the 77 component 

We now turn to quantifying the difference in {77) from figure [4] 
through the application of stellar population synthesis models. 
Since the change of 77 is not accompanied by a change in (, or the 
fractional change in high-£ galaxies, the shift is unlikely to be re- 
lated to recent star formation. Furthermore, this shift is relatively 
small, which implies composite models, such as an exponentially 
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decaying star formation history, may well blur any discriminating 
effect. Therefore, while the modelling of earl y-types can be push ed 
beyond relatively simple formation histories iRo eers et alj200^) . it 
is more robust in this case to consider simple stellar populations. 

The observations are compared to synthetic spectra through 
the equivalent widths of absorption line indices. The model 
spectra are extrac t ed fr om the s imple ste l lar po pulations of 
Bruzual & Chariot! d2003l) . using a IChabriej (2003) initial mass 
function. A detailed grid of models is constructed over a range 
of ages {1 ■ ■ • 14 Gyr} and metallicities —1.0 < \og(Z/Z@) < 
+0.35, at solar abundance ratios. We constrain the population pa- 
rameters with multiple age sensitive absorption features: H/3, H7, 
H<5, as well as the 4300 A absorption ba nd (G4300) and the metal- 
licity indicator [MgFe] dGonzalezll2003h . The absorption features 
are measured in each of the observed galaxy spectra, after being 
smoothed to the highest velocity dispersion within the respective 
bin. The eq uivalent widths (EW s) are estimated using the method 
outlined in Rogers et alj d2009h . in which the pseudo-continuum is 
determined by a boosted median value of the surrounding spec- 
trum. The specific absorption lines are targeted using a 20A win- 
dow. We use a standard maximum likelihood method to compare 
the observed and model absorption features. The model EWs are 
obtained in the same way, smoothing the spectra to the velocity 
dispersion considered. For each galaxy we use the probability- 
weighted age and metallicity, where the probability is defined as 
P(t, Z) ~ exp~° 5Ax , which gives a robust estimate of the SSP 
values (e.g. Gal lazzi et al.ll2005r) . The best fits give a reduced aver- 
age of (Xr) = 1.3. 

The average age and metallicity from each of the halo mass 
and velocity dispersion bins is used to investigate the difference 
seen in the average 77 value in an equivalent analysis to that shown 
in figure|4] However, we find no significant differences between the 
predicted SSP ages or metallicities with respect to host halo mass 
across all velocity dispersions. This illustrates the power of PCA to 
identify small differences in the spectra of elliptical galaxies, which 
might be difficult for model-based methods to identify, especially 
at low S/N. Therefore in order to assess the (second order) effects 
of environment on 77, we consider small perturbations with respect 
to the (first order) correlation with velocity dispersion. 

We define a fiducial relation between velocity dispersion and 
age/metallicity through the average parameters of each of the com- 
plete velocity dispersion bins i.e. including all halo masses. The 
SSP parameters of this relation are then offset with respect to the 
average age for each bin in velocity dispersion. The effect on the 
77 component of this change (computed directly on the models) is 
estimated and compared to the observed values - defining for each 
galaxy a Arj as the difference between the r\ value of the galaxy 
and that of the model corresponding to the same velocity disper- 
sion. The model values of 77 are derived in the same way as for the 
observed values, i.e. via projection onto the principal components 
followed by a rotation of the projected values. This rather simple 
analysis enables us to quantify the change in 77. Notice that we as- 
sume only a perturbation in age , as this is by far the most dominant 
effect reported in the lit erature dBernardi et alj|2003l ; lThomas et al.1 
l2005tlNelanetalj2005b . 

While we are only interested in identifying the magnitude 
of the shift of the lowest mass halo, as a check we compare 
our relationship of metallicity and age with velocity dispersion 
to those found by other authors. We quantify this relationship 
through the slope of a linear fit to the fiducial relation, which 
was found to be: A log(Z/Z©)/A log(cr) = 0.68 and A log(Age 
/Gyr) /A log (a) = 0.38 respectively. These values are com- 
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Figure 7. Comparison of the distribution of the central and satellite galaxy 
populations at the same velocity dispersion, with respect to average age 
((77); bottom) and to recent star formation (f(£ 0.5); top). The short 
dashed line is the (77) for the lowest mass halo in the whole sample. The 
thin solid line in the bottom panel is the reference line from fig[4] 



parable to those rep orted in iThomas et ail 1 20051) (0.55, 0.24), 



IClemens et afl d200fj) (0.76.~) and lNelan et alT il2005)(0.53.0.59). 

The a vs. A77 relationship is shown in figure|6]for the fiducial 
model obtained from the best fit SSP ages and metallicities (la- 
belled "Average"). For reference, we also show the relation when 
the age is shifted by 1 or 2 Gyr as labelled. We include in the fig- 
ure the observed values for the lowest halo mass (dashed line) and 
for the rest of the sample (thick solid line). The figure shows that 
early-type galaxies in the lowest density regions are on average 
about 1 Gyr younger than those in denser environments, a result 
that is consistent for a wide range of velocity dispersion. This re- 
sult indicates that galaxies residing in all but the lowest mass halos 
have similar stellar populations suggesting they are formed in sim- 
ilar ways at similar redshifts. This may be due to a true invariance 
across average/high density environments, such that above a certain 
density the formation process becomes uniform. 



3.4 Centrals and Satellites 

The star formation histories of galaxies sitting at the centre of the 
dark matter halos (centrals) and galaxies which orbit around the 
centre (satellites) are expected to depend in a different way with re- 
spect to the mass of their host halo. Specifically, satellites are more 
likely to be affected by the transformation mechani sms operating 
in the halo environment dvan den Bosch et al . 2008a b). Therefore, 
we analyse separately the effects of environment on centrals and 
satellites and compare the relative differences. 

The comparison of satellite and central early-type galaxies 
at the same velocity dispersion is shown in figure [7J The figure 
reveals that central galaxies (solid lines) have higher values of 
(rj) (bottom) and f(£ > 0.5) (top), indicative of younger aver- 
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age ages I Pasq uali et al . 2009b). and significant recent star forma- 
tion, respectively, when compared to satellit es of the same mass 
(long d ashe d lines). This is c onsist ent with Ivan den Bosch et all 
(2008b) and IWeinmann et al.l ( 120091) . who found that in terms of 
optical colours, centrals are bluer than satellite galaxies, indicating 
younger ages. Here we can see this extends to the fraction of ellip- 
tical galaxies which harbour small amounts of recent star formation 
as well. However, notice that this trend includes all halo masses. In 
§3.4.1 we show the effect of separating this central/satellite classi- 
fication with respect to halo mass. 

We have also plotted in the bottom panel of figure|7]the aver- 
age value of 77 for galaxies in the lowest mass halos (short dashed 
line). At the lowest halo masses all galaxies are central. Hence, 
this comparison allows us to rank the importance of host halo mass 
against the central/satellite nature. The figure shows that of these 
two properties, the mass of the halo dominates, i.e. centrals in low 
mass halos are younger than centrals in general. We focus on this 
aspect in the next two sub-sections. 

3.4.1 Effect on average age: r\ 

The average values of 77 are plotted for the two subgroups in fig- 
ure [8] It is evident from the figure that the analysis is complicated 
by the segregation imposed naturally by the halo mass on the cen- 
tral/satellite nature of the galaxy. The sample lacks central galaxies 
in the most massive halos (a limit imposed by the survey volume). 
On the other hand, at a given velocity dispersion, satellites cannot 
be found at low halo masses (a trivial constraint imposed by the 
halo mass). This means a comparison of the two populations can 
only occur within certain ranges of halo mass. Therefore the mid- 
dle panel of figure[8]shows the (rj) values only over a range of halo 
masses in which a significant number of both centrals and satel- 
lite galaxies exist. The top and bottom panels show the central and 
satellite galaxies, respectively, over their full halo mass range. 

The central galaxies (top), show an increase in the average 
value of {rj} at the lowest halo masses. This is identical to the trend 
shown in figure|4](notice centrals are the only type in the first halo 
mass bin). The result for more massive halos is consistent with the 
main sample. This suggests again that the effect of halo mass on the 
stellar populations is limited to masses below M H ~ 10 12 M e . On 
the other hand, satellite galaxies (bottom) show an increase in (77} 
in the group halo mass range 3 x 10 12 ^M H ^ 10 13 M Q , which is 
the lowest halo mass bin at which significant numbers of satellite 
galaxies exist. The comparison in the middle panel reveals that at 
the overlapping halo mass i.e. M H ~ (3 x 10 12 • ■ • 1O 13 )M , cen- 
trals have much lower values of (77). Thus the increase of (77) found 
in the population of satellite galaxies would have been hidden by 
the dominance of centrals in this mass bin. 

Given that higher 77 values are consistent with younger aver- 
age ages suggests that satellites are in fact younger, when compared 
to central galaxies of t he same velocity disp ersion and occupying 
halos of similar mass jPasauali et aT1l2009bl) . A KS test confirms 
that for this halo mass range, the second and third bins in a of the 
centrals and satellites are drawn from different distributions at a 
confidence level of Js98% and ^99%, respectively. While this is 
seemingly in contrast with the previous results, it is important to 
realise what is being compared in figure[8] A pure central - satellite 
split (as in figure |7]( will generally result in a comparison of cen- 
tral galaxies in low richness environments with satellites in higher 
density regions. This is the motivation of the split in many papers, 
under the assumption that the centrals of today are a reasonable ap- 
proximation to the progenitors of the current satellite population, 
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Figure 8. {77) values for the central and satellite galaxy populations within 
each bin in velocity dispersion mid halo mass. The central panel compares 
central and satellite galaxies over a range of halo masses for which sufficient 
numbers of both types exist. The halo mass range is as follows - given in 
log(M H /M ): grey dashed(11.80,12.15); grey dotted (12.15,12.50); grey 
solid(12.50,13.00); black dashed (13.00,13.50); black dotted(13.50,14.00); 
black solid(14. 00, 15.00). The thin solid line shown in all panels is the ref- 
erence line from fig|4] 



(e.g. §4. 1 Ivan den Bosch et ID |2008b)). Here we are comparing 
centrals in groups of significant size, to satellites in similar groups. 
We find a trend in three out of the four bins considered, although it 
is only found to be statistically significant in two. The effect is not 
visible at the higher velocity dispersions. 

Hence, we find that central galaxies 'evolve' faster than a sim- 
ilar galaxy that is accreted onto a group. This probably relates to 
the epochs over which environmental effects are efficient. Central 
galaxies within more massive halos will have interacted with other 
galaxies during the formation of the group and will also have ex- 
isted inside a more massive halo earlier than a galaxy of the same 
mass which is accreted int o a group halo . This result is consis- 
tent with the recent work of Ber nardi et alj (|2009j), who found that 
early-type brighest cluster galaxies were ~ lGyr older than the 
surrounding satellite early-type galaxies. 

Figure [8] reveals that the satellite galaxies in group halos of 
mass M H ~ 3 x 10 12 ■ ■ ■ 1O 13 M have values of (77) consis- 
tent with those of the central galaxies from halos of mass Mh ~ 
6 x 10 11 ■ ■ • 1O 12 M . This trend is consistent with the hierarchical 
build up of structures expected within the standard ACDM cosmol- 
ogy- 



3.4.2 Effect on recent star formation: £ 

Figure [9] shows the fraction of galaxies with values of £ consistent 
with recent star formation, with respect to central/satellite nature, 
in a similar manner to figure [8] The trends of the individual popu- 
lations of centrals (top) and satellites (bottom) are consistent with 
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Figure 9. Comparison of the distribution of central and satellite populations 
with respect to the fraction of galaxies with a f value above 0.5, which is 
consistent with recent star formation. The fraction is conditional i.e. calcu- 
lated for each bin separately.. The binning with respect to halo mass is the 
same as in figure [8] 

those found in the sample as a whole (see figure O. However a 
comparison of the fraction of high £ values between both centrals 
and satellites (middle) shows a significant difference. The central 
galaxy data (triangles) are consistently positioned at higher frac- 
tions of recent star formation even within the same halo mass range. 
This might not be particularly surprising: while satellites generally 
have their accretion of material stopped when they enter the halo, 
centrals do not, and so can s till accrete gas . This m ech anism is con- 
sistent with the scenario of iKavirai et alj d2009l) and lRogers et al.l 
(2009), whereby the recent star formation seen in elliptical galax- 
ies is fuelled by the accretion of small clouds of gas or satellites, 

3.5 Effects with Halo-Centric Radius 

Since there is evidence to suggest that the younger low halo mass 
galaxies are accreted as satellites, we might expect to see younger 
galaxies on the outskirts of groups. There is also the possibility that 
the halo-centric radius modulates the efficiency of environmental 
mechanisms. We investigate this aspect by looking at the relation- 
ship of rj and ( as a function of the projected distance from the 
luminosity weighted centre of the group. 

Shown in the right hand panel of figure [10] are the satel- 
lite galaxies split into subsamples according to velocity disper- 
sion and halo centric radius, (see fig. |2j, against the average val- 
ues of r\ for each bin. Since we are looking at galaxies across 
all group halo masses, the projected distance is scaled using the 
virial radius of the occupied halo. As found by pr evious work (e.g. 
Ivan den Bosch et al] |2008a; Wei nmann et alj|2009l) there is no sig- 
nificant trend. There is some indication that a reduction in the scat- 
ter (bottom panel) of the r\ values is seen for decreasing radii, al- 
though the effect is obviously far from robust. We note that this 



analysis is complicated both due to the degenerate nature of the 
projected distance as well as to the fact that smaller halos tend to 
be more concentrated. These effects are likely to blur the subtle 
relationships associated with elliptical galaxies and environment. 

We also analyse in the left hand panel of figure [10] the condi- 
tional fractions of galaxies which have £ ^ 0.5 as a function of the 
scaled projected distance from the luminosity weighted centre. An 
inspection of the figure indicates a general trend of increasing re- 
cent star formation fractions with projected radius, consistent with 
the previous comparison of centrals and satellites. The stripping of 
a satellite surrounding gas may well explain both the trend seen 
here and in figure [9] Note that the recent star formation fractions 
seen at the outskirts of the groups are higher than those for the gen- 
eral satellite sampl e at the same stellar mass. This is consisten t with 
the results seen in Ferreras et ai] d2006h : iRogers et all &wfo ; and 
lRogersetalJd2009h , where we show that galaxies in medium den- 
sity environments show increased fractions of recent star formation 
possibly due to interactions. The large error bars are a function of 
the small numbers which exist away from the group centres, yet the 
trend is consistent across all mass bins. 

To conclude, we find that the effects with halo centric distance 
are seen mainly in terms of £ i.e. younger subpopulations. As a re- 
sult of being accreted into the halo an elliptical galaxy will have any 
residual star formation halted. On the other hand, we find that the 
average age of the bulk of the stellar component does not depend 
on the distan ce from the centre of the halo . This is to be expected if, 
as argued by , van den Bosch et al. (2008a), a large fraction of early- 
type galaxies transition onto the red sequence as centrals. It may 
be that the m ain effects expecte d are better analysed within cluster 
samples (e.g. iNelan et al.ll2005h . 



4 DISCUSSION AND CONCLUSIONS 

Using the large ^7,000 strong sample of early-type galaxies de- 
scribed in Roge rs et al.l d2007h we have investigated the effect of 
environment, measured through the mass of the dark matter halo 
that hosts each galaxy. We make use of the previously derived PCA 
rotated projections, 77 and £ , sensitive to the average properties and 
recent star formation, respectively, to identify small differences in 
their stellar populations. We find that the star formation histories 
are mainly a function of velocity dispersion. This is shown in fig- 
ure|3]where a considerable difference is evident between the effects 
of velocity dispersion and group halo mass. When we separate out 
the sample to remove the mass-environment degeneracy, the stellar 
populations across most of the halos are indistinguishable. There- 
fore, while the majority of this paper has focused on the small dif- 
ferences observed in the elliptical galaxy population, the main con- 
clusion is that the effect of environment on the stellar populations 
of early-type galaxies is extremely limited. However, we do find 
that small but nonetheless interesting differences can be detected 
with respect to environment. 

First of all, galaxies within halos with the lowest masses are 
estimated to be ~ 1 Gyr younger than those in the rest of the sam- 
ple. This offset is consistent with the fact that galaxies in low 
mass halos are exclusively centrals. Thus, they represent the lower 
end of the environmental density scale and so our result compares 
favourably with the general v iew that galaxies in low density re- 
gions are 1—2 Gvr younger {Bernardi et al.lll998l; iThomas et al.l 
l2005l : lBernardi et alJl200rj ; ISanchez-Blazquez et alj2006l) . 

We use in this paper a catalogue of groups to estimate the 
dark matter halo masses, allowing us to compare more naturally 
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Figure 10. The r\ and f components for the satellite subsample are shown as a function of distance from the luminosity weighted centre of the galaxy group, 
where the distance is scaled using R v ; r . (Left:) The main panel shows the average r\ per bin in R pro j/R v i r , where adjoined points correspond to the same range in 
velocity dispersion (labelled in km/s). The smaller panel below shows the standard deviation within each bin. (Right:) Fraction of galaxies with a £ value above 
0.5 (i.e. those with signatures of recent star formation). The fraction is computed within each bin separately. The error bars are the fractional Poisson errors. 



with theoretical modelling results. On e of the more recent discov- 
eries is that many observations (e.g. iBinnevI 20041. ICroton et al.l 
120051 : iDekel & Birnboinjl2006l : ICattaneo et alj|2008h can be well 
explained by assuming a critical halo mass above which the supply 
of cold gas is stopped by virial shock heating teirnboim & Dekell 
l2003l:lKeres et al.l2005l ; lDekel & Birnb oim 2006). The critical halo 
mass i s found to be Mh ~ 10 12 M Q . Modelling bv lCattaneo et al.l 



(2008) show that galax ies above and bel ow this value are consistent 
with the analysis o f lThomas et al] fc005) with respect to galaxies in 
high/low density regions (see Cattaneo et al. figure 6). It is possible 
that we see this divide more explicitly here, since the lowest halo 
bin of this sample is the only one to contain galaxies in halos below 
or close to the critical mass halo considered here. Our sample gives 
a consistent ~ 1 Gyr age difference between galaxies in halos with 
masses above and below this critical mass. We do not find a grad- 
ual trend from the lowest mass halo upwards, a result that may be 
caused by a combination of effects. The li mited timescale of evo- 
lutionary signatures on the optical spectra lHarke retal1l2006l) and 
the importance of the individual galaxy mass and the nature of the 
group build up or the low quality of the data, may mean that the 
difference in SFH is only visible in the lowest halo mass where the 
effect is strongest. 

The average age is not the only parameter that varies across 
the sample. The fraction of galaxies with a high value of £ (con- 
sistent with recent star formation) is higher in the four lowest mass 
halos. The reason for the drop above Mh ~ 3 x 10 13 /i _1 Mq is 
not entirely obvious. However, we note that this is the same halo 
mass at which a decline i s observed in the op tical AGN popu- 
lation jGilmour et al]|2007l ; IPasauali et al1 r2009). possibly giving 
way to radio-mode AGN. We also note that gas stripping processes 
are more effective in more massive halos. The modelling results of 
Simh a et al. | j2009h indicate that satellites of medium and low mass 
i.e. those containing the majority of the recent star formation, have 
their accretion stopped efficiently at halo masses M H ~ 1O 14 M . 



We also observe a decrease of recent star formation in the satel- 
lite population as a whole, as well as with decreasing halo centric 
radius, suggesting that such process works across all halo masses 
although with differing efficiency. 

We also investigated the differing effects of environment on 
the population of central and satellite galaxies. These galaxies are 
dominant in low and high mass halos, respectively, which implies 
a limited overlap on the parameter space spanned by Mh and a. 
In the range of halo mass over which a comparison is possible, 
the two populations were found to have similar stellar populations 
at high velocity dispersion (a ^200 km/s). However, at lower 
values (150 km/s ^ a ^ 200 km/s), (77) is higher for satellite 
galaxies, as expected for younger ages. Hence, the stellar popu- 
lations of central galaxies were formed earlier. Environmental ef- 
fects preferentially act on satellites whereas central galaxies mainly 
move onto the red sequence as the res ult of a merger or when it s 
halo mass surpasses the critical mass o f IDekel & Birnboim! (2006 ) . 
Central galaxies generally sit at the p eaks of the dark matter den- 
sity distribution (Berli nd et al. 2003). In halos with mass Mh ~ 
3 x 10 13 /i _1 M Q , centrals are within groups containing significant 
numbers of additional (satellite) galaxies. Therefore, they are likely 
to have been subject to a higher level of interactions at earlier times 
and possibly had gas heate d by infalling material and satellites 
dKhochfar & Ostr iker 2008). Furthermore, if we assume that the 
central galaxy is the core member of the grou p, we would expect 
these galaxies to have crossed the threshold of IDekel & Birnboim! 
(2006) earlier than a similar satellite, which would have been ac- 
creted later. This would imply that at least some satellites will have 
become ellipticals after being accreted onto the group, creating the 
lower mean ages. Therefore, we have the emerging picture that 
while central galaxies are quenched on average at earlier times they 
retain or accrete small amounts of gas with which to form small 
amounts of stars. 
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